Consistent description of nuclear charge radii and electric monopole transitions 
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A systematic study of energy spectra throughout the rare-earth region (even-even nuclei from 
58 Ce to 74 W) is carried out in the framework of the interacting boson model (IBM), leading to an 
accurate description of the spherical-to-deformed shape transition in the different isotopic chains. 
The resulting IBM Hamiltonians are then used for the calculation of nuclear charge radii (including 
isotope and isomer shifts) and electric monopole transitions with consistent operators for the two 
observables. The main conclusion of this study is that an IBM description of charge radii and 
electric monopole transitions is possible for most of the nuclei considered but that it breaks down 
in the tungsten isotopes. It is suggested that this failure is related to hexadecapole deformation. 
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I. INTRODUCTION 

The structural properties of excited 0"*" states in de- 
formed even-even nuclei have been the subject of a long 
controversy. According to the geometric model of Bohr 
and Mottelson [l|, a nucleus with an elipsoidal equi- 
librium shape may undergo oscillations of two different 
types, l3 and 7. The first type of vibration preserves axial 
symmetry while the second allows excursions toward tri- 
axial shapes. These vibrations should be combined with 
rotations exhibited by the deformed, vibrating nucleus 
to yield a rotation-vibration spectrum. The geometric 
model of Bohr and Mottelson, therefore, predicts an axi- 
ally deformed nucleus to display a spectrum of rotational 
bands built on top of vibrational excitations. The low- 
est in energy is the ground-state rotational band with 
K^ = 0"*" (i.e., of which the projection of the angular 
momentum on the axis of symmetry is zero) which cor- 
responds to no intrinsic excitation. Next in energy are 
the /3- and 7-vibrational bands which correspond to one 
intrinsic excitation of the /? or 7 type, characterized by 
a rotational band with K"^ — 0+ or 2+, respectively. At 
higher energies still, the geometric model predicts rota- 
tional bands built on multiple excitations of (3 and/or 7 
phonons. 

While 7-vibrational bands are an acknowledged feature 
of deformed nuclei, such is not the case for /3- vibrational 
bands. Confusion arises because excited 0^ states in nu- 
clei can be of many different characters, such as pairing 
isomers Q , two-quasi-particle excitations (S] , or so-called 
intruder states that arise through the mechanism of shape 
coexistence [J] . A careful analysis of the observed prop- 
erties of excited O"*" states seems to indicate that very 
few indeed satisfy all criteria proper to a /3-vibrational 
state [3- In particular, although this observation is ob- 
fuscated by a lack of reliable data, very few O"*" states 
decay to the ground state by way of an electric monopole 
transition of sizable strength 6] , as should be the case for 
a /? vibration 7]. It is therefore not surprising that al- 
ternative interpretations of excited 0"*" states in deformed 
nuclei, either as pairing isomers (see, e.g., Refs. @,@]) or 



through shape coexistence and configuration mixing (see, 
[lO|) have gained advocates over recent years. 
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The purpose of this paper is to examine to what ex- 
tent a purely collective interpretation of nuclear 0"^ lev- 
els is capable of yielding a coherent and consistent de- 
scription of observed charge radii and electric monopole 
transitions [Il|. As noted above, very few measured 
electric monopole transitions satisfy the criteria proper 
to a matrix element from the ground-state to the /3- 
vibrational band and the present attempt therefore might 
seem doomed to failure. However, collective excitations 
of nuclei can also be described with the interacting bo- 
son model (IBM) of Arima and lachello |12h14| |. where 
they are modeled in terms of a constant number of s 
and d (and sometimes g) bosons which can be thought of 
as correlated pairs of nucleons occupying valence shell- 
model orbits coupled to angular momentum £ = and 
2 (and 4), respectively. One of the advantages of the 
IBM is that a connection with the shell model [l^| as 
well as with the geometric model |16l - [l8J | has been estab- 
lished. In particular, one of its dynamical symmetries, 
the SU(3) limit [13], displays energies reminiscent of the 
rotation-vibration spectrum of the geometric model. It 
has also been shown, however, that the first-excited 0+ 
state in the SU(3) limit of the IBM has not exactly a 
/3-vibrational character but is a complicated mixture of 
intrinsic /3 and 7 vibrations [19j. The main purpose of 
this paper is to show that a collective interpretation of 
excited 0"*" states with the IBM is not inconsistent with 
the electric monopole data, as observed in the rare-earth 
region. 

The outline of this paper as follows. In Sect. |TT] the 
ground is prepared by discussing charge radii and electric 
monopole transitions in the context of different models 
that are applied to even-even nuclei in the rare-earth re- 
gion from Ce to W in Sect. lIIII A qualitative explanation 
of the failure of this approach in the W isotopes is offered 
in Sect. IIVI bv invoking effects of hexadecapole deforma- 
tion. Finally, Sect. IVl summarizes the conclusions of this 
work. 



II. CHARGE RADII AND ELECTRIC 
MONOPOLE TRANSITIONS 

Electric monopole (EO) transitions between nuclear 
levels proceed mainly by internal conversion with no 
transfer of angular momentum to the ejected electron. 
If the energy of the transition is greater than 2mcC^ 
(where mo is the mass of the electron), they can occur 
via electron-positron pair creation. A less probable de- 
excitation mode which can proceed via an EO transition is 
two-photon emission. It is not a priori clear why a con- 
nection exists between charge radii and EO transitions. 
In fact, the argument is rather convoluted and we begin 
this section by recalling it. The argument can be general- 
ized to effective operators, leading to a relation between 
charge radii and EO transitions which forms the basis of 
the present study. 
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where the first summation runs over all A nucleons while 
the second and third summations run over neutrons (v) 
and protons (tt) only, and where Cn (ep) is the effec- 
tive neutron (proton) charge. If bare nucleon charges 
are taken (en — and ep — e), the summation extends 
over protons only and the original expression ^ is re- 
covered. If equal nucleon charges are taken (en = ep), 
Eq. ([4]) is appropriate for the matter radius. 

In the approximation cr « 0, the starting expres- 
sions ^ and (121) for the nuclear EO transition strength 
and the mean-square charge radius are identical (up to 
the constants R^ and Z). It is therefore natural to fol- 
low the same argument as used for the charge radius for 
the construction of a generalized EO transition operator, 
leading to the expression [2l| 



A. Relation between effective operators for charge 
radii and electric monopole transitions 

The total probability for an EO transition between ini- 
tial and final states |i) and |f) can be written as the prod- 
uct of an electronic factor fl and a nuclear factor p(EO), 
the latter being equal to [20| 
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with R = tqA^'^ and where the summation runs over the 
Z protons (tt) in the nucleus. The coefficient a depends 
on the assumed nuclear charge distribution but in any 
reasonable case it is smaller than 0.1. The second term 
in Eq. ^ therefore can be neglected if the leading term 
is not too small [20j. In this approximation we have 
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On the other hand, the mean-square charge radius of 

a state |s) is given by 
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This is an appropriate expression insofar that a realis- 
tic A-body wave function is used for the state |s). For 
the heavy nuclei considered here the construction of such 
realistic wave function is difficult and recourse to an ef- 
fective charge radius operator T{r'^) should be taken. In 
particular, if neutrons are assigned an effective charge, 
the polarization of the protons due to the neutrons is 
'effectively' taken into account, giving rise to changes in 
the charge radius (r^) with neutron number. The gen- 
eralization of the expression ^ can therefore be written 
as 
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In terms of this operator, the dimensionless quantity 
p(EO), defined in Eq. ([2]) and referred to as the monopole 
strength, is given by 



p(EO) = 



(f|r(EO)|i 
ei?2 



(6) 



Since the matrix element ^ is known up to a sign only, 
usually p2(E0) is quoted. 

The basic hypothesis of the present study is to assume 
that the effective nucleon charges in the charge radius 
and EO transition operators are the same. If this is so, 
comparison of Eqs. (jlj and ([5]) leads to the relation 



f (EO) = (eniV - 



epZ)f(r2). 
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This is a general relation between the effective operators 
used for the calculation of charge radii and EO transitions, 
which is applied throughout this study. 



B. Charge radii and electric monopole transitions 
in the interacting boson model 

Equation ([7]) can, in principle, be tested in the frame- 
work of any model. This endeavor is difficult in the con- 
text of the nuclear shell model because realistic wave 
functions, appropriate for the calculation of (r^) or EO 
matrix elements, are hard to come by for heavy nuclei. 
As an alternative we propose here to test the implied 
correlation with the use of a simpler model, namely the 
IBM [13 - [l4| . This requires that all states involved [i.e., 
|s) in Eq. ^, and |i) and |f) in Eq. ^] are collective in 
character and can be described by the IBM. 

In the IBM-1, where no distinction is made between 
neutron and proton bosons, the charge radius operator is 



taken as the most general scalar expression that i 
in the generators of U(6) J22l] . 
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where Ni^ is the total boson number (the customary no- 
tation N is not used here to avoid confusion with the 
neutron number), hd is the d-boson number operator, 
and a, 77 are parameters with units of length^. The first 
term in Eq. (|8]), (r^)c, is the charge radius of the core nu- 
cleus. The second term accounts for the (locally linear) 
increase in the charge radius due to the addition of two 
nucleons {i.e., neutrons since isotope shifts are considered 
in this study). The boson number Ny^ is the number of 
pairs of valence particles or holes (whichever is smaller) 
counted from the nearest closed shells for neutrons and 
protons. If the bosons are particle-like, the addition of 
two nucleons corresponds to an increase of A^b by one 
and a is positive. In contrast, if the bosons are hole-like, 
the addition of two nucleons corresponds to a decrease 
of A^b by one and a is negative. Therefore, care should 
be taken to change the sign of a at mid-shell \22\. The 
third term in Eq. (|S]) stands for the contribution to the 
charge radius due to deformation. It is identical to the 
one given in Ref. [22] but for the factor 1/N\j. This fac- 
tor is included here because it is the fraction {hd)/N\^ 
which is a measure of the quadrupole deformation (/S^ in 
the geometric collective model) rather than the matrix 
element {fid) itself. Since the inclusion of 1/A^b is non- 
standard, also results without this factor will be given in 
the following, that is, with the charge radius operator 



T'(r2) = (r2)e + a'iVb + ?7'"d- 



(9) 



Once the form of the charge radius operator is deter- 
mined, that of the EO transition operator follows from 
Eq. ([T]). In the IBM-1 the EO transition operators are 
therefore 



The occurrence of the absolute value \a\ is due to the 
interpretation of the bosons as pairs of particles or holes, 
as discussed above. Isomer shifts are a measure of the 
difference in charge radius between an excited (here the 
21) state and the ground state, and are given by 
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Similar formulas hold for the charge radius operator (|9]) 
in terms of the parameters a' and 77' . For ease of notation, 
the superscript {A) in the isotope and isomer shifts shall 
be suppressed in the following. 



C. Estimate of the coefficients a and r] 

Although the coefficients a and rj in Eq. ([8]) will be 
treated as parameters and fitted to data on charge radii 
and EO transitions, it is important to have an estimate of 
their order of magnitude. The term in a increases with 
particle number and therefore can be associated with the 
'standard' isotope shift. This standard contribution to 
the charge radius is given by [23| 
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The term in 77 stands for the contribution to the nuclear 
radius due to deformation. For a quadrupole deformation 
it is estimated to be |23|1 
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where P2 is the quadrupole deformation parameter of the 
geometric model. 

The estimate of lal follows from 



T(EO) = (enAT + epZ)/; 
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Since for EO transitions the initial and final states are 
different, neither the constant (r^)c nor aA^b in Eq. (|S]) 
or a'A^b in Eq. © contribute to the transition, so they 
can be omitted from the EO operators. 

Two other quantities can be derived from charge radii, 
namely isotope and isomer shifts. The former measures 
the difference in charge radius of neighboring isotopes. 
For the difference between even-even isotopes one finds 
from Eq. ^ 
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which for the nuclei considered here {A ~ 150) gives |a| w 
0.2 fm^. The estimate of ry can be obtained by associating 
the deformation contribution (J15l) with the expectation 
value of (n-d)o+ in the IBM. This leads to the relation 
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where use has been made of the approximate correspon- 
dence /32 ~ {4,Nh/^A)^/^T^2 between the quadrupole de- 
formations P2 and 82 in the geometric model and in the 
IBM, respectively [2J|. The relation dTT]) yields the esti- 
mate 
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For typical values of TVb ^ 10 and A ^ 150 this gives 
a range of possible rj values between 0.25 and 0.75 fm^, 
corresponding to weakly deformed (/32 ^1) and strongly 
deformed (/32 ~ V^) nuclei, respectively. 

Similar estimates can be derived for the coefficients in 
the alternative form ^ of the charge radius operator. 
The estimate for \a'\ is identical to Eq. ((T6)) while the 
one for 77' differs by a factor N^,. 



D. Estimate of the effective charges 

The consistent definition of operators for charge radii 
and EO transitions leads to the introduction of neutron 
and proton effective charges for both operators, as op- 
posed to previous treatments where this was done for EO 
transitions only. This opens the possibility to obtain an 
estimate of the effective charges en and Cp from data for 
charge radii which are widely available. 

A simple estimate can be obtained in a Hartree-Fock 
approximation with harmonic-oscillator single-particle 
wave functions. For the Hartree-Fock ground state |gs) 
of a nucleus, simple counting arguments of the degen- 
eracies of the three-dimensional harmonic oscillator (see 
Sect. 2.2 of Ref. [25[), lead to the following expression for 
the expectation value of the neutron part of the charge 
radius operator: 
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where b is the oscillator length of the harmonic oscilla- 
tor for which the estimate b = iJ I^Z-^I'^^-^I'^tqA^I^ is 
used [25|] . Equation (fT9|) and its equivalent for the proton 
part of the charge radius operator, can then be combined 
to yield 

^'^^«^" 5 "^^ e„7V-HepZ ' ^^"^ 

The separate effective charges Cn and Cp cannot be ob- 
tained from a fit of this expression to the data on charge 
radii, but the ratio en/sp can. There exists, unfortu- 
nately, a strong dependence of this ratio on the value 
used for Tq, as will be discussed in Subsect. IIIICl 



Assume that the initial and final levels in the EO tran- 
sition are orthogonal mixtures of some states \oi\J) and 
\oiiJ), so that 



|iJ) = a./laiJ) + 6j|a2J), 
|fj) = aj\a2J) -bj\a\J). 



(21) 



The labels a\ and a2 refer to a different intrinsic struc- 
ture for the two sets of states (typically two bands) , the 
members of which are additionally characterized by their 
angular momentum J . The aj and bj are mixing coeffi- 
cients that satisfy Oj -h 6j = 1. If one assumes that no 
EO transition is allowed between states with a different 
intrinsic structure, (ai J|T(E0)|a2 J) ~ 0, it follows that 



(fJ|T(EO)|iJ) 
« a.ibj Ua2J\f{m)\oL2J) - (ai J|T(EO)|ai J) ' 
= aj6j(e„iV-t-epZ)((r2)„,j-(r2)„^j), (22) 

where the last equality is due to Eq. ((T]). Furthermore, 
it must be hoped that the difference in (r^) appearing in 
Eq. ([^ can be identified with a measured isotope shift 
A(r^) between nuclei somewhere in the neighborhood 
and, therefore, that the ground states of two nuclei in 
the neighborhood can be identified with the unmixed in- 
trinsic structures a\ and a-i- Finally, it must be assumed 
that this difference in {r^) does not depend significantly 
on J or, equivalently, that isomer shifts are identical for 
the intrinsic structures a\ and a^- With this first set of 
assumptions the following relation holds: 
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which reduces to the result of Wood et al. Q] if bare nu- 
cleon charges are taken. According to this equation the J 
dependence of the EO strength is contained in the coeffi- 
cients a J and 6j, which can be obtained from a two-state 
mixing calculation. An additional input into the calcu- 
lation, therefore, is the relative position in energy of the 
intrinsic structures ai and a2 and the size of the mixing 
matrix element (itself assumed to be independent of J). 
These are the essential ingredients of the calculations re- 
ported by Kulp et al. [10| of which the results are quoted 
below. 



E. Relating charge radii and electric monopole 
transitions through configuration mixing 

The idea of correlating observed EO transitions with 
charge radii is not new and, in particular, was already 
proposed by Wood et al. [0]. These authors tested this 
proposal in the context of a model assuming mixing be- 
tween two different coexisting configurations. Since in 
Sect, mil results are quoted for the few nuclei where this 
model has been applied, a brief reminder of the method 
is given in this subsection. 



III. SYSTEMATIC STUDY OF NUCLEI IN THE 
RARE-EARTH REGION 

To test the relation between charge radii and EO tran- 
sitions, proposed in the previous section, a systematic 
study of all even-even isotopic chains from Ce to W is car- 
ried out. This analysis requires the knowledge of struc- 
tural information concerning the ground states and ex- 
cited states which here is obtained by adjusting an IBM-1 
Hamiltonian to observed spectra in the rare-earth region. 



A. Hamiltonian and energy spectra 

All isotope series in the rare-earth region from Z — 58 
to Z = 74 have the particularity to vary from spher- 
ical to deformed shapes and to display systematically 
a shape phase transition. Such nuclear behavior can 
be parametrized in terms of a simplified IBM-1 Hamil- 
tonian which can be represented on the so-called Cas- 
ten triangle [2g, as demonstrated for rare-earth nuclei 
by Mc Cutchan et a l. ul} . Alternatively, as shown by 
Garcia- Ramos et al. [2al , the same region of the nuclear 
chart can be described with the full IBM Hamiltonian. 
The latter approach is adopted here and a general one- 
and two-body Hamiltonian is considered which is written 
in multipolc form as |22| 



H = ednd+aoP+-P-+aiL-L+a2Q-Q+a3T3-T3+a4T4-T4, 

„ ^ (24) 
where hd is the d-boson number operator, P-|_ (P-) is 
a boson-pair creation (annihilation) operator, L is the 
angular momentum operator, and Q, T^, and T4 are 
quadrupole, octupole, and hexadecapole operators, re- 
spectively. (For the definitions of these operators in 
terms of s and d bosons, see Sect. 1.4.7 of lachello and 
Arima |22i].) If only excitation and no absolute ener- 
gies are considered, the expression ([M)) defines the most 
general Hamiltonian with one- and two-body interactions 
between the bosons in terms of six parameters. The pa- 
rameter X in the quadrupole operator Q can be arbitrar- 
ily chosen (as long as it is different from zero) and is fixed 
here to — -\/7/2 for all nuclei. 

Although reasonable results are obtained with constant 
parameters for an entire chain of isotopes, the spherical- 
to-deformed transition is better described if at least one 
parameter is allowed to vary with boson number TVb. 
Garcia-Ramos et al. Q followed the procedure to vary 
the d-boson energy Cd with iVb. A different method is 
followed here by allowing the variation of the coefficient 
a2 associated with the quadrupole term. To reduce the 
number of parameters in the fit, we assume that 02 de- 
pends linearly on the quantity N,yNTr/{N^ + A^^) where 
N„ (Nt^) is half the number of valence neutrons (pro- 
tons) particles or holes, whichever is smaller. The ar- 
gument for introducing such a dependence is related to 
the importance of the neutron-proton interaction as the 
deformation-driving mechanism in nuclei [29|. The pa- 
rameter 02 can then be decomposed into two terms as 
follows: 
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The parameters in the Hamiltonian ([M)) are deter- 
mined from a least-squares fit to levels of the ground- 
state band with K'^ = 0+ and those of two more bands 
with K'' = 0+ and A'" = 2+ ('quasi-/3' and 'quasi-7' 
bands). Since two-quasi-particle excitations do not be- 
long to the model space of the IBM-1, states beyond 
the backbend cannot be described in this version of the 



model and for this reason only levels up to J'^ = 10"'" are 
included in the fit. Similarly, near closed shells, excita- 
tions might be of single-particle character and, therefore, 
nuclei with A^ < 84 are excluded from the energy fit. 
Nevertheless, the Hamiltonian ((24)) with the parametriza- 
tion P5|) allows the extrapolation toward the A^ = 84 
isotopes, which is needed for some of the isotope shifts 
calculated in the following. The total number of nuclei 
included in the fit is 78 and the total number of excited 
levels is 846. The parameters are summarized in Ta- 
ble U as well as the root-mean-square (rms) deviation a 
for each isotope chain which typically is of the order of 
100 keV. 

Figure [T] illustrates, with the examples of samarium, 
gadolinium, and dysprosium, the typical evolution of the 
energy spectrum of a spherical to that of a deformed 
nucleus which is observed for every isotope series studied 
here. In each nucleus are shown levels of the ground-state 
band up to angular momentum J^ — 10"*" as well as the 
first few states of the excited bands, together with their 
experimental counterparts, if known. 

For a given nucleus the choice of the subset of col- 
lective states that should be included in the fit is often 
far from obvious. While members of the ground-state 
and the 7-vibrational bands in a deformed nucleus are 
readily identified, this is not necessarily so for the '/3- 
vibrational' band. As a result, guided by the EO tran- 
sitions discussed in the next subsection, the calculated 
first-excited K^ — 0"*" band is associated in some nuclei 
not with the lowest observed K'^ = 0"*" band but with a 
higher- lying one. This is the case for 168-172YI3 a,nd also 
for i^^Er where the fourth J'^ = 0+ level at 1934 keV 
has been recently identified as the band head of the (3- 
vibrational band on the basis of its large EO matrix ele- 
ment to the ground state (SOJ • 

Although the agreement with the experiment can be 
called satisfactory, it should be noted from Table U that 
it is obtained with rather large parameter fluctuations 
between the different isotopic chains. This presumably is 
so because the parameters are highly correlated and small 
changes in the fitted data give rise to large fiuctuations in 
some of the parameters. The main purpose of this calcu- 
lation, however, is not to establish some parameter sys- 
tematics with the Hamiltonian ([M]) but rather to arrive 
at a reasonably realistic description of the spherical-to- 
deformed transition. This will enable a simultaneous cal- 
culation of charge radii and EO transitions, as discussed 
in the next two subsections. 



B. Isotope and isomer shifts 



Isotope shifts A(r^), according to Eq. p^. depend 
on the parameters \a\ and 77 in the IBM-1 operator ([H]). 
These parameters are expected to vary smoothly with 
mass number A, according to Eqs. pB]) and (fTSl) . The 
estimate p6p neglects, however, the microscopic make-up 
of the pair of neutrons which varies considerably from Ce 



TABLE I: Parameters in the Hamiltonian (|24|l and the rms deviation a, in units of keV. 
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FIG. 1: Experimental (points) and calculated (lines) energy levels in the Sm, Gd, and Dy isotopes: (a) ground-state band up 
to r = 10+, (b) first-excited K" = 0+ band up to J" ^ 4+ , and (c) first-excited K" = 2+ band up to r = 4 + . 



TABLE IL The parameters |a| and \a'\ in the charge radius 
operators ([8]) and ((9|, in units of fm^, for the different isotope 
series. 



Isotope 


Ce 


Nd 


Sm 


Gd 


Dy 


Er 


Yb 


Hf 


W 


\a\ 
\a'\ 


0.22 
0.23 


0.24 
0.25 


0.26 
0.26 


0.13 
0.09 


0.15 
0.12 


0.15 
0.12 


0.11 
0.09 


0.10 
0.11 


0.11 
0.15 



to W. It is therefore necessary to adjust \a\ for each iso- 
tope series separately, and the resulting values are given 
in the first row of Table |IT1 The value of 7/, on the other 
hand, is kept constant for all isotopes, 77 = 0.50 fm^. The 
parameters thus derived are broadly consistent with the 
estimates obtained in Subsect. Ill CI 

Similar remarks hold for the alternative parametriza- 



Nd 




^y 



Gd 



Dy 



Er 



o o O o 



0.25 


0.25 


0.25 



0.25 

_ 

t^ 0.25 

- 

0.25 


0.25 


0.25 



0.25 



90 100 110 

Neutron number A'^ 

FIG. 2: Experimental (points) and calculated (lines) isotope 
shifts A(r^), in units of fm^, for isotopic chains in the rare- 
earth region from Ce to W. The full (dashed) lines are for 
the charge radius operator ((81 [(O], with parameters give n in 
Table lU and in the text. Data are taken from Ref. [sj for 
Ce, from Ref. [H for Nd, Sm, Dy, Er, and Ybjrom Ref. [H 
for Gd, from Ref. [H for Hf, and from Ref. ^ for W. 



tion ^ of the charge radius operator. The values of 
I a' I are given in the second row of Table HIl while rj' = 
0.05 fm^. 

The resulting isotope shifts are shown in Fig. [5] Only 
mariginal differences exist between the two sets of calcu- 
lations, with and without the factor l/A'b in the charge 




radius operator. The main reason for preferring the 
form ([S]) is that it lacks a kink in A(r^) at mid-shell. Al- 
though this seems to occur in the Hf data, it is unlikely 
that the observed kink is associated with a maximum 
of the boson number at mid-shell. Given the similar- 
ity between the two sets of calculations, the subsequent 
comments are valid for both. 

The peaks in the isotope shifts are well reproduced in 
all isotopic chains with the exception of Yb. The largest 

peaks occur for 152-150gj^^ 154-152q^^ ^^^ 156-154j3y^ 

that is, for the difference in radii between A^ = 90 and 
A^ = 88 isotopes. The peak is smaller below Z = 62 
for Ce and Nd, and fades away above Z = 66 for Er, 
Yb, Hf, and W. The calculated isotope shifts broadly 
agree with these observed features but there are differ- 
ences though. Notably, the calculated peak in the Sm 
isotopes is much broader than the observed one, indi- 
cating that the spherical-to-deformed transition occurs 
faster in reality than it does in the IBM-1 calculation. 
Also, it would be of interest to determine the character 
of the transition in the Nd isotopes: the present IBM-1 
calculation predicts it to be rather smooth but data in the 
deformed region of the transition are lacking to confirm 
this behavior. Likewise, the IBM-1 calculation features 
a very fast transition for the Gd isotopes with a sharp 
peak at A^ = 88 but isotope-shift data are lacking for the 
nuclei in the spherical region of the transition. 

The feature of peaking isotope shifts is related to the 
onset of deformation which is particularly sudden (as a 
function of neutron number) for the Sm, Gd, and Dy iso- 
topes. This can be quantitatively understood as a result 
of the subshell closure at Z = 64 [30|, combined with 
the strongly attractive interaction between neutrons in 
the lj^ft.9/2 orbit and protons in the lT:hiii2 orbit. If the 
occupancies of the neutrons in the 1 1^/19/2 and of the pro- 
tons in the l7r/iii/2 orbit are both low, as is the case for 
A^ < 92 and Z < 64, the nucleus is expected to be spher- 
ical. As soon as one of the two orbits becomes signif- 
icantly occupied, the strong neutron-proton interaction 
will induce occupancy of the partner orbit and an onset 
of deformation. 

A further test of the calculated charge radii is ob- 
tained from isomer shifts 5{r'^), depending only on rj [see 
Eq. 113])] or rj' . The isomer shifts that are known experi- 
mentally are listed in Table IIIII The data are more than 
30 years old and often discrepant. Nevertheless, a clear 
conclusion can be drawn from the isomer shifts measured 
in the Sm and Gd isotopes: they are easily an order of 
magnitude smaller in the deformed than they are in the 
spherical region. In spite of the extreme sensitivity of 
this effect, a quantitative description is obtained of the 
isomer shifts in the Gd isotopes. For the Sm isotopes only 
a qualitative agreement is found since the experimentally 
observed drop in isomer shift between ^'^^Sm and ^^^Sm 
is stronger than what is calculated in the IBM-1. This in- 
dicates that the spherical-to-deformed transition is faster 
in reality than it is in the calculation, in line with what 
can be concluded from the isotope shifts. 
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TABLE III: Experimental and calculated isomer shifts 5{r'^), 
in units of 10" fm , in the rare-earth region. 







5(r2> (10-^ 


fm^) 






Isotope 


Thl" 


Th2'' Th3" 


Expt 


Ref 


i™Sm 


80 


72 


49.6 


2.6 


m 


is^Sm 


37 


37 19 


25. 
19. 


7. 


[38] 
[39] 








14. 


1. 


[40] 








12. 




m 








18. 


4. 


[37] 


is^Sm 


6 


7 


1.1 


0.8 


[42] 


is^Gd 


28 


31 


15. 


2. 


[40] 








19. 


6. 


[43] 








18.5 


2.5 


[44] 








20.0 


3.9 


[45] 


i56Gd 


4 


5 


2.6 


0.8 


[46] 








0.1 


1.7 


[44] 








4.3 


3.7 


[45] 


i5»Gd 


3 


4 


0.4 


0.3 


[47] 








1.5 


0.8 


[48] 








-4.0 


1.2 


[44] 








0.3 


3.3 


[45] 


i^^Gd 


1 


1 


0.3 


0.8 


[48] 








-1.4 


3.1 


[44] 








-2.1 


3.2 


[45] 


170Yb 


4 


6 


1.2 


0.3 


[46] 








1.7 


0.6 


[49] 








1.20 


0.40 


[50] 


172 Yb 


3 


5 


0.41 


0.20 


[50] 


174Yb 


2 


4 


1.0 


0.45 


[49] 








-0.44 


0.19 


[50] 


176Yb 


2 


4 


-0.17 


0.10 


[50] 


182W 


5 


6 


6.0 
-0.6 
-0.2 




[51J 
[41J 
[52] 


184W 


8 


9 


0.16 




[52] 








0.5 


0.3 


[53] 


186W 


17 


18 


0.14 




[52] 



"With the char; 
''With the char; 
'^Configuration 



;ge radius operator JSj. 
-ge radius operator iQl l. 
mixing calculation of Ref. [lCl| . 



From the preceding analysis the following picture 
emerges. All considered isotopic chains exhibit an evolu- 
tion from a spherical to a deformed shape which, at the 
phase-transitional point, is characterized by a peak in the 
isotope shifts. The height of the peak is proportional to 
the suddenness of the transition. This effect is a direct 
consequence of the increase in the mean-square radius of 
a nucleus due to its deformation. The IBM-1 is able to 
provide an adequate description of this transitional be- 
havior. By adjusting the charge radius operator of the 



IBM-1 to the observed height of the peak in the isotope 
shifts, a first estimate of the parameter rj (or rj') is ob- 
tained. Its value follows more directly from isomer shifts 
since only one parameter enters this quantity but, unfor- 
tunately, data are scarce and often unreliable. The choice 
rj = 0.50 fm^ (or rj' = 0.05 fm^) is a compromise between 
the value obtained from a fit to A(r^) of all isotopes and 
the one from S{r^) in the Gd isotopes. The question is 
now whether this value of ry (or 77') reproduces the EO 
transitions observed in the rare-earth nuclei. 



C. Electric monopole transitions 

The calculation of the matrix elements of the EO tran- 
sition operator pUj) or (|lip requires the knowledge of the 
effective charges e^ and ep. In principle, an estimate of 
the ratio Cn/ep can be obtained by fitting the expres- 
sion (|^D|) to the available data on charge radii in the 
rare-earth region. The minimum in the rms deviation 
is shallow though and, furthermore, the correlation be- 
tween ro and Cn/cp is strong. In other words, a slightly 
different choice of tq gives an almost equally good fit 
to the charge radii of 58 < Z < 74 nuclei but with a 
significantly different ratio en/cp. A reasonable choice 
of parameters, close to the optimum set, corresponds to 
tq = 1.24 fm, Bn — 0.50e, and Cp = e. 

In Table IIVI the available EO data in the rare-earth re- 
gion are compared with the results of this calculation. 
The two choices of EO transition operator, Eqs. (fTO)) 
and (llip . again yield comparable results. An overall 
comment is that the present approach succeeds in re- 
producing the correct order of magnitude for p^(EO), in 
particular in the Sm, Gd, and Dy isotopes. However, 
some discrepancies can be observed in heavier nuclei and 
especially concern ^^^Yb and i82-i84^^ ^ possible ex- 
planation is that the /9^(E0) measured for these nuclei is 
not associated with collective states. This seems to be 
the case in ^^^Yb where several p^(EO) have been mea- 
sured none of which is large. Only in the W isotopes does 
it seem certain that the observed EO strength is consis- 
tently an order of magnitude smaller than the calculated 
value. It is known that these nuclei are in a region of 
hexadecapole deformation [55{ and this may offer a qual- 
itative explanation of the suppression of the EO strength, 
as argued in the next section. 

While in a spherical vibrator there is no appreciable EO 
strength from the ground state to any excited 0+ state, 
this is different in a deformed nucleus which should ex- 
hibit large p'^(EO)s from the ground-state towards the 
/3- vibrational band [l|,|3l- As a consequence, one predicts 
an increase in the EO strength as the phase-transitional 
point is crossed. This seems to be confirmed in the few 
isotopic chains where data are available. Adopting a 
simple, schematic Hamiltonian, von Brentano et al. [56| 
showed that also in the IBM-1 sizable EO strength should 
be observed in all deformed nuclei. The present IBM-1 
calculation is in qualitative agreement with this geomet- 



TABLE IV: Experimental and calculated p^(EO) values in the 
rare-earth region. 





Transition 


J 




p^(EO) X 10 


i 


Isotope 


Thl" 


Th2' Th3" 


Expf^ 


isoSm 


740 ^ 





7 


6 


18 2 




1046 -^ 334 


2 


16 


13 


100 40 


^^^Sm 


685 -> 





52 


52 72 


51 5 




811 -> 122 


2 


41 


41 77 


69 6 




1023 -^ 366 


4 


29 


29 84 


88 14 




1083 -^ 





2 


2 


0.7 0.4 




1083 -> 685 





47 


47 


22 9 


is^Sm 


1099 -^ 





41 


49 


96 42 


is^Gd 


615 -^ 





68 


68 


63 14 




931 -^ 344 


2 


77 


77 


35 3 


i54Gd 


681 ^ 





84 


102 


89 17 




815 -5> 123 


2 


66 


80 


74 9 




1061 ^ 361 


4 


38 


46 


70 7 


is^Gd 


1049 -^ 





44 


64 


42 20 




1129 ^ 89 


2 


41 


59 


55 5 


ISSQd 


1452 -> 





30 


51 


35 12 




1517 -^ 79 


2 


27 


45 


17 3 


158Dy 


1086 -^ 99 


2 


42 


70 


27 i2 


leoDy 


1350 ^ 87 


2 


28 


56 


17 4 


le^Er 


1171 -^ 102 


2 


38 


64 


630 460 


164E, 


1484 -s> 91 


2 


24 


48 


90 50 


««Er 


1460 -^ 





9 


20 


127 60 


170Yb 


1229 ^ 





32 


72 


27 5 


172Yb 


1405 -> 





30 


76 


0.2 0.03 


174jjf 


900 -^ 91 


2 


32 


71 


27 W 


176Hf 


1227 -> 89 


2 


15 


38 


52 9 


178Hf 


1496 ^ 93 


2 


32 


72 


14 3 


182W 


1257 ^ 100 


2 


45 


77 


3.5 0.3 


184w 


1121 -^ 111 


2 


52 


75 


2.6 0.5 



"With the EO transition operator (IIOI I. 
'With the EO transition operator ((TT}. 
'^Configuration-mixing calculation of Ref. [lO|. 
"^From Ref. [5411 for J = 0, except ^^^Sin and ^^^Er which are from 
Ref. [33 ; from Ref. Q for J ^ 0. 



ric picture and with the resuhs of von Brentano et al. 
Nevertheless, it should be pointed out that, systemati- 
cally, the calculated p^(EO; Oj — >■ 0^) in Table HV] dimin- 
ishes once the phase-transitional point is crossed. In the 
Gd nuclei, at least, this behavior seems to be borne out 
by the data. It indicates that the first-excited 0+ state in 
the IBM-1 has not simple a /3- vibrational character but 
has a more complicated structure [3, [I^l • 

Tables Hill and HV] also show the results of a 



configuration-mixing calculation for ^^^Sm [l^]. This ap- 
proach leads to a quantitative, detailed description of the 
data. Results of similar good quality are obtained for 
^^''Gd [57|. However, in view of the employed method- 
ology, as explained in Subsect. Ill El it seems difficult to 



make systematic calculations of EO properties of nuclei 
with this model. 



IV. EFFECT OF g BOSONS ON ELECTRIC 
MONOPOLE TRANSITIONS 



An obvious extension of the sd-IBM is to include a cor- 
related pair of higher angular momentum for which the 
most natural choice is the g boson with i = 4. Many arti- 
cles have been published over the years where the role of 
the g boson has been investigated in detail, for which we 
refer the reader to the review of Devi and Kota [58| . The 
sdg-lBM has been used in the interpretation of struc- 
tural properties of nuclei in the rare-earth region. For 
example, properties of the is^-ieoQ^^ isotopes, including 
energy spectra and E2, E4, and EO transitions were inter- 
preted in the framework of the sdg-IBM-1 [50|. The con- 
clusion of this particular study, namely that the g boson 
is indispensable for the explanation of the character of 
some K'^ — 4~^ bands, was confirmed in transfer-reaction 
studies, see, e.g., Burke et al. [60,|6i|. Other examples 
of studies of rare-earth nuclei in the sdg-IBM-l include 
i68Er [13, "6-i58sin [63(, and 144-150^^ ^^. 

In this section the possible influence of hexadecapole 
deformation or, equivalent ly, of g bosons on EO tran- 
sitions is discussed. The spherical-to-deformed shape 
phase transition in the sdg-IBM-1 corresponds to a tran- 
sition between the two limits U(5)(g)U(9) and SU(3) ^ . 
The following schematic Hamiltonian is adopted: 



H = edfid + egfi 



a"'g 



(26) 



where Q^ is the SUsrfg(3) quadrupole operator [67 



)^ = [st X J+dt xS](2) 



11 

14 



2Mtxd](2) 



(27) 



+ y[rf^ xg + g^ X 



J](^)-^V55[5t^^l(^^ 



x.9j^ 



For a convenient description of the phase transition, an- 
other parametrization of the Hamiltonian (j26p can be 
introduced in terms of A and Q (sometimes referred to as 
control parameters) which are related to e^, e^, and k by 



fid 



c 



4A^b(l-C)' 



(28) 



where A'b is now the total number of s, d, and g bosons. 
The Hamiltonian (1^51) then becomes 



H^c 



(i-C)(«rf + An,)-^g.g 



(29) 



where c is a scaling factor. The U(5) (8) U(9) limit is 
obtained for ^ = whereas the SU(3) limit corresponds 
to C = 1. By varying C, from to 1 one will cross the 
critical point ^c ~ 0.5 at which the spherical-to-deformed 
transition occurs. 
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In the sdg-IBM-l the charge radius operator is 



T(r2) = (r2), + a^b + r?^+7^, (30) 

while the EO transition operator is 



T(EO) = (e„iV + epZ) ( r/^ 



'iVb 






(31) 



which are straightforward extensions of the expres- 
sions ([8]) and (fTO| . Again, the total number operator 
-^b — TT-s '^ fi-d + fig docs uot Contribute to the EO transi- 



tion and is not included in the operator (1311) . 

Analytic expressions can be derived for the matrix el- 
ements of the operators fi^, n^, and fig for the limiting 
values of C in the Hamiltonian (P^ . They are known for 
arbitrary angular momentum J but for simplicity's sake 
results are quoted for J = only. In the U(5) (g) U(9) 
limit they are trivial, 



{Ot\ns\0', 



N^Sa: 



(0+|n,|0+) = 0, 
(0^l"9|0+> = 0. 



(32) 



In the SU(3) limit one finds for the ground-state expec- 
tation values. 



(0+|n,|0+) 
{Ot\nd\Ot) 
(0+|n,|0+) 



JVb(4A^b + 1) 
5(4A^b-3) ' 

16(iVb - l)A^b(4A^b + 1) 
7(4iVb-3)(47Vb-l) ' 
64(iVb - l)A^b(2iVb - 3) 
35(4iVb-3)(4iVb-l) ' 



(33) 



and for the transition matrix elements from the ground 
to the first-excited 0+ state, 



(0+|n.|0+) 
{Ot\nd\0+) 
(0+|nJO+) 



2(iVb-l)iVb(2iVb-l)(4iVb + l) 



3(4iVb-5)2(4iVb-3) 
2(iVb-l)iVb(2Afb-l)(4iVb 



1/2 



13)2(47Vb + l) 



1/2 



3(4Arb - 5)2(4iVb - 3)(4iVb - 1)^ 

96 r 2(JVb - l)iVb(2A^b - 3)2(2JVb - l)(4iVb + 1) 
35 [ 3{ANy, - 5)2(4Arb - 3)(4iVb - 1)2 

I 



1/2 



(34) 



TABLE V: Matrix elements in the classical limit of SU(3) in 
the sd-IBM-1 and the sdp-IBM-1. 





(0+|n,|0+> 


(0+|n,|0+) 


IBM 


i=G 1=2 e=A 


e 


= 


£ = 2 


e = 4: 


sd 
sdg 


iiVb In, - 

14 8 

-iV. -N, -N, 


2 

3\ 
2 
5\ 


In, 

' 2 

In, 

' 3 


2 In, 

3V 2 

2 In, 

7V 3 " 


24 In, 

35 V 3 



It is instructive to compare these results to the cor- 
responding ones in the sd-IBM-1 which is done in Ta- 
ble |V] in the classical limit iVb — > 00. Considering first 
the expectation values of n£ in the ground state, one 
notes that d bosons are dominant in OJ" both in the sd- 
and sdg-IBM-l and that the contribution of g bosons 
in the sdg-IBM-l is fairly modest. One therefore does 
not expect a significant effect of the g boson on the nu- 
clear radius, and this should be even more so away from 
the SU(3) limit for a realistic choice of boson energies, 



< e^ < Cg. In other words, the schematic Hamilto- 
nian (|26p captures the obvious feature that effects of de- 
formation on the nuclear radius are mainly of quadrupole 
character, and that hexadecapole deformation plays only 
a marginal role. In terms of model calculations it also 
means that the parameter 7 in the operator (1301) is iff 
determined from radii because the expectation value of 
fig in the ground state is small. Nevertheless, on physi- 
cal grounds one expects rj as well as 7 to be positive and 
of the same order since both the quadrupole and hex- 
adecapole deformation have the effect of increasing the 
nuclear radius. 

Turning to the (Oj^ln^jOj ) matrix elements, one notes 
first of all that in the sd-IBM-1 only i = 2 contributes 
to the 0^ -^ 0^ EO transition with a coefficient 77 fixed 
from the isotope shifts. As discussed in the preceding 
section, this typically leads to a large p2(E0) from the '/3- 
vibrational' to the ground band. The situation is, how- 
ever, drastically different in the sdg-IBM-l. It is seen 
from Table |V] that the matrix elements of fig is larger 
than that of fid and of different sign. Therefore, while 
changes in the nuclear radius due to the g boson are ex- 
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<oii«rf|or>.d 
<orinrfior>sd, 

<0||n,|0|>sd, 




the sd^-IBM-l equal to t] in the sd-IBM-l all results of 
the preceding section concerning radii are reproduced. 

In the sd-IBM-l as well as in the sd^-IBM-l a sharp 
increase in (Aj'^lndjOj)^ is observed around Cc ~ 0.5, see 
Fig. Mjo)- Up to that point, ( < 0.5, there is essen- 
tially no contribution to p^(EO; 0^ — > 0^) from the g bo- 
son. Consequently, all sd-IBM-l results up to the phase- 
transitional point are not significantly modified by the g 
boson. As can be seen from Fig. El^b), in the deformed 
regime this is no longer true since, in the sd^-IBM-l, a 
sharp decrease of (0]*"|nrf|Oj)^ occurs at C ~ 0.6 and, fur- 
thermore, (0^|ng|Oj)^ rapidly increases beyond ( « 0.5 
and dominates (Of |nd|Oj)^ for ( > 0.7. The explanation 
of this dominance is that (Of |nd|Oj) changes sign before 
reaching its value in the SUsdg(3) limit, in agreement 
with the analytical results quoted above. 



4- <0||«rf|0j>2, 
(0||«rf|0j>2 



(b) 



;dg 



3- 



(0||«g|Oj>?j. 




FIG. 3: The matrix elements (a) {Of\he\0^) and (b) 

(0]^|n^|Oj)'^ for £ = 2 and ^ = 4 in the spherical-to-deformed 
transition of the sd-IBM-1 and the sdg-lBM-1. In the 
sd-IBM-1 the transition is from U(5) to SU(3) and in the 
sdg-mM-1 from U(5) (g) U(9) to SU(3) with A = eg/ea = 1.5. 
The number of bosons is A^b ~ 8. 



pected to be negligible, one cannot rule out its significant 
impact on p^(EO; Oj — >■ Of) in deformed nuclei. 

This argument can be made more quantitative by 
studying the spherical-to-dcformed shape transition of 
the Hamiltonian (1261) . The matrix elements of hd and fig 
can be calculated for arbitrary C with the numerical code 
ArbModel [6g]. A reasonable choice for the ratio of boson 
energies is A = 1.5. The results are shown in Fig. |3] and 
compared to the matrix elements of hd calculated for the 
U(5)-to-SU(3) transition in the sd-IBM-I. Figure M.a.) 
confirms the dominance of the d boson in the ground state 
of deformed nuclei both in the sd- and sdg-lBM-1. More- 
over, the expectation value of fid varies with ^ in very 
much the same way in both models. In fact, for the en- 
tire transition the relation {0'l\fis\0'^)sd ~ (Of |ns|Of )sdg 
approximately holds, meaning that by choosing 77 -I- 7 in 



V. CONCLUSIONS 

In this paper we proposed a consistent description of 
nuclear charge radii and electric monopole transitions. 
The ingredients at the basis of such a description are (i) 
the derivation of a relation between the effective opera- 
tors describing nuclear charge radii and electric monopole 
transitions, (ii) the mapping of these operators from the 
shell model to the interacting boson model, (iii) the de- 
scription of spectroscopic properties of chains of isotopes 
through the shape-transitional point with the interact- 
ing boson model, and (iv) the assumption that initial 
and final states in the considered electric monopole tran- 
sitions have a collective character and can be adequately 
described with the interacting boson model. 

The validity of this approach was tested with an ap- 
plication in even-even nuclei in the rare-earth region 
(58 < Z < 74) which systematically display a spherical- 
to-deformed transition. This transitional behavior could 
be successfully reproduced with the interacting boson 
model and was shown to be correlated with peaks in the 
isotope shifts, as observed at the phase-transitional point. 
In particular, the correlation between the suddenness of 
the shape transition and the height of the peak in the 
isotope shift could be correctly reproduced by the model. 
With the charge radius operator determined in this way 
from isotope and isomer shifts, an essentially parameter- 
free and systematic calculation of electric monopole tran- 
sitions in the rare-earth region could be undertaken. The 
observed electric monopole strengths were reproduced to 
within a factor 3, except in the isolated case of ^^^Yb 
and in the W isotopes. As a possible explanation for the 
failure of the approach in the latter isotopes, the role of 
hexadecapole deformation or, equivalently, of the g bo- 
son was explored in a schematic model. It was concluded 
that the effect of the g boson is marginal on charge radii 
but can be strong on electric monopole transitions. 

We are aware that our explanation of electric monopole 
strength is based on a geometric picture of the nucleus, in 
contrast to an alternative explanation in terms of shape 
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coexistence and configuration mixing. While there are 
undoubtedly regions of the nuclear chart {e.g., Sr, Zr, and 
Mo isotopes) where the latter mechanism is needed to 
explain the observed electric monopole strength, we have 
presented here a comprehensive analysis of this quantity 
in rare-earth nuclei that lends support to the geometric 
interpretation. 

Many thanks are due to Stefan Heinze who helped us 
with the numerical calculation involving the g boson. We 
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